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Gravitino dark matter from Q-ball decays 
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Affleck-Dine baryogenesis, accompanied by the formation and subsequent decay of Q-balls, can 
generate both the baryon asymmetry of the universe and dark matter in the form of gravitinos. 
The gravitinos from Q-ball decay dominate over the thermally produced population if the reheat 
temperature Tr < 10 7 GeV. We show that a gravitino with mass ~ 1 GeV is consistent with all 
observational bounds and can explain the baryon-to-dark-matter ratio in the gauge-mediated models 
of supersymmetry breaking for a wide range of cosmological and Q-ball parameters. Moreover, 
decaying Q-balls can be the dominant production mechanism for m 3 / 2 < 1 GeV gravitinos if the 
Q-balls are formed from a (B — L) =0 condensate, which produces no net baryon asymmetry. 
Gravitinos with masses in the range 50 eV < m 3 / 2 < 100 keV produced in this way can act as warm 
dark matter and can have observable imprint on the small-scale structure. 

PACS numbers: 12.60.Jv, 95.35.+d 



I. INTRODUCTION 

Supersymmetric extensions of the standard model pre- 
dict the existence of baryonic and leptonic Q-balls [1], 
which may be stable, or may decay into fermions. In 
the early universe, large Q-balls are abundantly produced 
from the fragmentation of a flat direction condensate. In 
case of gaugc-mcdiatcd supersymmetry breaking, these 
Q-balls can be stable [2] and exist today as dark mat- 
ter [3]. The cosmology and astrophysical implications of 
Q-balls have been studied by a number of authors [4-40] . 
Here we study the cosmological implications of small- 
charge, unstable Q-balls, which may not exist at present, 
but whose formation and decay in the early universe 
could produce the population of dark matter particles 
consistent with observations. 

A priori, one could expect the baryon-to-dark matter 
(BDM) ratio, ilb/^DM to be very different from one if 
dark matter and baryonic matter arise from very different 
processes. The observational fact that ftb/^DM = 0.17, 
is within one order of magnitude from unity is intriguing. 
This may be an indication that baryons and dark matter 
share a common physical origin [41]. Most mechanisms 
of dark matter production are completely unrelated to 
baryogenesis and therefore implicitly assume that this 
ratio is an accident. In the case of Affleck-Dine cosmol- 
ogy [23, 42], dark matter can originate from the same 
process as baryogenesis, which suggests a possible expla- 
nation to the fact that Qb and Qdm are not drastically 
different. Models along these lines have been constructed 
for both gravity-mediated and gauge-mediated super- 
symmetry breaking models. In gauge- mediated models 
the partial evaporation of baryons off the Q-ball are the 
source of baryons, while the remaining Q-balls act as 
dark matter [4, 5, 43]. In gravity-mediated models the 
late decay of Q-balls produces baryons and non-thermal 
LSP dark matter [6, 7]. 

We will show that decays of unstable Q-balls in the 
gauge-mediated scenario can serve as the dominant and 
non-thermal source of dark matter in the form of graviti- 
nos, hence relating dark matter and baryon densities in 



a new way. In our model a ~ 1 GeV gravitino is favored. 
In conventional scenarios, the gravitinos are produced in 
a combination of thermal processes and NLSP decays, 
^3/2 ~ 1 GeV. Such a population of gravitinos can be 
dark matter, but the scenario requires a finely tuned re- 
heating temperature, 10 6 GeV < T R < 10 7 GeV [44]. In 
models of gaugc-mcdiatcd SUSY breaking in which the 
flat direction dynamics are properly taken into account, 
the reheat temperature is often significantly below this 
desired value [45, 46]. Our production mechanism ex- 
tends the viability of gravitino dark matter to low re- 
heating temperature. 

While a somewhat similar model of gravitino pro- 
duction from Q-ball decays has been studied before in 
the context of gravity-mediated SUSY breaking [47] , the 
GeV scale of gravitino mass is theoretically disfavored in 
gravity-mediated scenarios. Moreover, the model consid- 
ered in Ref. [47] required a high reheating temperature 
and neglected some important thermal processes in cal- 
culating the Q-ball decay temperature. 

Enqvist and McDonald [6, 7] have discovered a very 
simple way of creating the correct amounts of dark mat- 
ter and baryon number simultaneously. The crucial ob- 
servation is that the squarks q in unstable Q-balls decay 
into quarks q and the LSP via q — > LSP + q. For every 
baryon number lost by the Q-ball at least three LSPs are 
created, Nlsp > 3. Since Q-ball formation takes place 
non-adiabatically the fraction of baryon number trapped 
in the Q-balls, /s, is generally maximal /b ~ 1 [48]. 
Then one finds 
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(1) 



Since Q-balls are a large conglomeration of squarks, it 
should not be surprising that in an i?-parity preserv- 
ing theory their decay necessarily creates comparable 
amounts of LSPs and baryons. In the original model, En- 
qvist and McDonald [6, 7] considered gravity-mediation, 
in which the LSP dark matter was the neutralino. How- 
ever a GeV neutralino is theoretically disfavored and now 
experimentally ruled out by LEP. To maintain the pre- 
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diction of eq. (1), Roszkowski and Seto [49] extended this 
idea to the case of axino LSP in gravity mediation which 
naturally accounts for the GeV scale. However the model 
is most natural in the gauge-mediated case, which al- 
ready has a light LSP, namely, the gravitino. This is 
intriguing because gravitinos with 0(1 GeV) masses are 
not only natural in gauge-mediated supersymmetry mod- 
els, but this scale is particularly attractive in its ability 
to solve several cosmological and phenomenological prob- 
lems [50, 51]. Moreover, although large gauge-mediated 
Q-balls can be stable, they are unstable for Q < 10 12 
at zero temperature and evaporate by the present-day 
time from thermal effects for Q < 10 18 [43]. It is these 
Q-balls wc focus on as the source of the correct BDM 
ratio. Moreover, decaying Q-balls can be the dominant 
production mechanism for m 3 / 2 < 1 GeV gravitinos if the 
Q-balls are formed from a (B — L) = condensate which 
produces no baryons. In this case of course the novel 
feature is a new mechanism of gravitino production, and 
the BDM ratio cannot be accounted for. 



II. REHEATING TEMPERATURE 

To avoid gravitino overclosure, the reheating tempera- 
ture must not exceed the following upper bound [52, 53]: 



m 3/2 \ ( TeV 
M^ 3 



(2) 



where Mq 3 is the gluino mass. For a GeV gravitino the 
upper bound on the reheating temperature is 10 s GeV. A 
more stringent constraint however comes from requiring 
that the late decays of the NLSP not disrupt successful 
big bang nucleosynthesis: Tr < 10 7 GeV [54, 55]. The 
reheat temperature in the case of stop or sneutrino NLSP 
is less constrained than in the case of stau NSLP: 
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III. 



GRAVITINO ABUNDANCE FROM Q-BALL 
DECAYS 
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Gravitinos freeze out at a very high temperature [44, 
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For the 1 GeV gravitinos with TeV-scale gluinos the 
freeze-out temperature is Tf ~ 2 x 10 13 GeV. As long as 
Q-balls decay well below this temperature, the produced 
gravitinos do not thermalize. Notice that we must have 
Tr -c Tf to avoid gravitino overclosure. However, grav- 
itinos are produced from scatterings in the hot plasma of 
MSSM particles. 



We now compute the Q-ball decay temperature and 
show that the produced gravitinos do not thermalize. At 
finite temperature the effective potential along a flat di- 
rection is [3, 4, 43] 



V(<f>,T)=m 4 (T) log 1 



m 2 (T) 



(5) 



where m(T) = max(M s ,T), Ms is the supersymmetry 
breaking scale, and T is the temperature. The mass of 
a squark inside the Q-ball is u> <~ m(T)/Q 1 / 4 . Notice 
then that while all Q-balls with Q < 10 12 are unstable at 
low-temperature (T < Ms ~ 1 TeV), larger Q-balls may 
decay as well via q — > G + q at high temperature. 

At finite temperature, both thermal evaporation of Q- 
balls and accretion of global charge onto Q-balls are pos- 
sible [4, 5, 8, 10, 15, 43]. The fate of Q-balls is, to some 
extent, determined by the difference between the chemi- 
cal potential of the surrounding plasma /U p i asma and inside 
the Q-ball /j,q ~ ui. For sufficiently small charge Q-balls 
at high temperature, the Q-ball chemical potential dom- 
inates [iq 3> Hpiasma, thus leading to the evaporation of 
such Q-balls [4], so that all Q-balls with charges Q < 10 18 
would have evaporated by the present era [43]. 

At high temperatures the effects of diffusion are impor- 
tant, while at lower temperatures thermal evaporation 
dominates. The temperature at which this transition oc- 
curs is T* i=a M s « 1 TeV for small charge Q-balls. Ka- 
suya and Kawasaki [43] have shown that at early times 
(T > Tr), the charge transfer rate is 
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(6) 



Thus the charge evaporated from a Q-ball during this 
epoch is 



AQ = 3MTj\ 
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where we have evolved the temperature from to Tr, 
assuming T > T R . Thus for Tr < 10 7 GeV Q-balls with 
charges Q > 10 11 survive until after reheating. 

Q-balls decaying before reheating decay at the temper- 
ature (T > Tr) 
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Since the Q-balls evaporated in this era are small Q < 
10 11 and the Tr < 10 7 GeV, the decay temperature 
Td <C Tf, which ensures that the gravitinos produced 
do not thermalize. 

The Q-balls in the range 10 11 < Q < 10 18 can survive 
until well after reheating. If they decay when diffusion 
is still dominant, then their decay temperature is (T* < 
T D < T R ) 
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The largest Q-ball that decays in this regime is Q ~ 10 16 . 
One can show that the Q-balls in the range 10 16 < Q < 
10 18 all decay very close to T D w 10 2 GeV. 

Thus the Q-balls in the range 10 12 < Q < 10 18 de- 
cay after reheating, producing baryons and non-thermal 
gravitinos. The gravitino abundance comes potentially 
form three sources: thermal production, non-thermal 
production from the NLSP decays, and from Q-ball decay 
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However in the case of charged slepton NLSP, their decay 
contributions to the gravitinos is subdominant compared 
to thermal processes for m 3 / 2 ^ 10 GeV [44]. We there- 
fore focus on the comparison between Y§ and Yg p . The 
bound on the reheating temperature already ensures that 
the first two terms will not lead to overclosure. The num- 
ber density owing to Q-ball decays is 



Y§ = N G f B 



n b 



(11) 



Using rib/s « 5 x 10 10 with Nq « 3 and /g w 1, we 
have Yq ~ 10~ 9 . The abundance from thermal processes 
is [44, 52] 
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Thus for a GeV gravitino the Q-ball production mecha- 
nism is dominant (Yq > Yq p ) as long as Tr < 10 7 GeV. 
This automatically guarantees that the gravitino and 
the NLSP do not affect the successful predictions of 
conventional cosmology. Gravitinos produced from ei- 
ther thermal processes or Q-ball decay extend the via- 
bility of gravitino dark matter to reheating temperatures 
Tr < 10 7 GeV, though only when Q-ball decay is dom- 
inant can the BDM ratio be successfully accounted for. 
Additionally the dependence on the specific flat direc- 
tion is weak, although we need (B — L) ^ to prevent 
sphalerons from washing-out any of the produced B and 
L asymmetry. 

It has been pointed out [56] that long-lived, negatively 
charged electroweak-scale particles X~ can overproduce 
Li 6 via the reaction ( 4 HeX- 1 ) + D -> Li 6 + X". Such 
overproduction of Li 6 constrains the lifetime of the X~ 
to be t < 5 x 10 3 s. This constraint is directly relevant 
to gauge-mediated models in which charged sleptons are 
often the NLSP. However we note that the NLSP lifetime 
is given by 
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where m T is the stau mass. For GeV-scale gravitinos 
this implies the very weak bound m T > 60 GeV. Thus 
no significant constraints from BBN are imposed on GeV 
gravitinos. 



IV. BEYOND GcV-SCALE GRAVITINOS 

Even if we abandon the desire to account for 
(1](,/J]dm), the gravitinos produced from Q-ball decay 
can dominate other production mechanisms for to 3 / 2 < 
GeV and Tr < 10 7 GeV. Since in this case we are not 
attempting to explain t]b ~ 10~ 10 , the Q-balls must 
be formed from a (B — L) = flat direction, since 
sphalerons will destroy any baryon asymmetry. The 
gravitino yield is still described by Eq. 11, but now 
rib is replaced by nm-L) which can be much larger 
V(b-l) > 10~ 10 [28, 40]. This feature allows the grav- 
itino yield from Q-ball decays to exceed the yield from 
other sources. Now however there is a lower-bound on 
the decay temperature of the Q-balls, To > T ew , so that 
the produced baryon/lepton asymmetries are erased. As 
we have seen, however, this constraint is not difficult to 
satisfy. 



V. FREE-STREAMING LENGTH 

Although the GeV-mass gravitinos from Q-ball de- 
cays constitute a cold dark matter candidate, smaller 
mass gravitinos can have a cosmologically relevant free- 
streaming length. One can estimate the free-streaming 
length by calculating the present-day free-streaming ve- 
locity v [44, 57] 
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The free-streaming length limits from the Lyman-a 
forrest constrain Xps ^ 0.59 Mpc [44, 58], which trans- 
lates to the limit v < 0.05 km/s on the free-streaming 
velocity. At the time of decay the momentum of the 
gravitino is roughly 
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which yields the free-streaming velocity 
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(16) 



where T is the present-day temperature, Tjj is the Q-ball 
decay temperature and m p is the proton mass. 

We note that since small Q translates into high de- 
cay temperature (early times), the small Q-balls have 
a chance of producing warm dark matter. The earlier 
produced gravitinos that are warmer. For the GeV- 
scale gravitino, the present-day free-steaming velocity is 
< 10~ 9 km/s, which corresponds to conventional cold 
dark matter. 
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FIG. 1: Here we show the parametric dependence of the free- 
streaming length Xfs on Q-ball charge Q and gravitino mass 
OT3/2- The blue region is excluded because it gives rise to hot 
dark matter with Xfs ^ 1 Mpc. 



The gravitinos with masses less than GeV may have a 
cosmologically relevant free-streaming length. In Fig. 1 
we plot the limits on the Q-ball produced gravitinos for 
Xfs — 1 Mpc, 100 kpc, 10 kpc, and 1 kpc. There are 
different limits on Xfs from different observations. There 
are limits based on the Lyman-a forrest at z = 3 (Xfs ~ 
0.59 Mpc) [58], while reionization has very different limits 
(Xfs ^ 0.94 Mpc) [59]. 

It has been argued in fact, that thermal gravitino dark 
matter is now ruled out for 7773/2 < 16 eV based on 
Lyman-a limits at z w 2 — 3 [60]. However this con- 
straint is only directly applicable to the thermally pro- 
duced gravitino dark matter. If instead the reheating 
temperature is small Tr < 10 7 GeV as we have argued 
here, then the results of [60] are evaded. A number of 
similar proposals for evading Lyman-a bounds exist. For 
example, entropy production from messenger decay [61] 
and the late decay of the relic neutralino population both 
produce gravitinos that are cool enough to evade Lyman- 
a constraints [51]. Both of these solutions however prefer 
a particular gravitino mass: 1 keV in the case of messen- 
ger entropy production [61]; and 7773/2 ~ (1 — 10) GeV 
from late neutralino decay [51]. As can be seen from 
Fig. 1 however, Q-ball produced gravitinos are viable over 
a much wider range of masses: 50 eV < m 3 / 2 ^3 100 keV. 

Moreover, observations of the motions of stars in dwarf 
spheroidal galaxies show evidence of a non-zero free- 
streaming length [62]. Observations therefore seem to 
prefer a free-streaming length close to the Lyman-a lim- 



its. Thus for a relatively wide range of parameters Q-ball 
decay can produce phcnomcnologically viable warm dark 
matter. 

VI. CONCLUSIONS 

As long as Tr < 10 7 GeV, Q-balls with baryon num- 
bers 10 12 — 10 18 decay after reheating and preserve the 
successful Enqvist-McDonald relation relating baryonic 
and dark matter energy densities 



^3/2^ \ _ / j 
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Here we took Ng = 3 and Jb ~ 1 , consistent with numer- 
ical simulations [48]. Thus the gravitino mass consistent 
with the correct BDM ratio is 



1.6 GeV 
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2.0 GeV 
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(18) 



It is instructive to compare the production of grav- 
itinos from Q-ball decays, the thermal freeze-out pro- 
cesses, and the non-thermal late decays of NLSPs. In 
particular, the three possibilities have different implica- 
tions for the reheat temperature. If the thermally pro- 
duced population dominates the relic abundance, then 
the reheating temperature must be finely-tuned to give 
^t-Dhih 2 ~ 0.11. If the contribution from NLSP decay 
dominates, then the reheating temperature must be rel- 
atively small Tr <C 10 9 GeV and generally requires a 
larger gravitino mass 7713/2 ^ 10 GeV. Finally, the pro- 
duction from Q-ball decays requires a small reheating 
temperature Tr ^ 10 7 GeV and small gravitino mass 
TO 3/2 ^ 1 GeV. In the special case with the gravitino 
mass 777,3/2 ~ 1 GeV one can account for the correct 
baryon-to-dark matter ratio. 

Finally, if the Q-balls form from a (B — L) = Affleck- 
Dine condenstate(c/. the models in Refs. [28, 40]), the 
decays of Q-balls produce no net baryon asymmetry that 
can survive the sphaleron transitions, but these decays 
can still act as a source of non-thermal gravitinos. In 
this case, if the gravitino mass is in the range 50 eV < 
m 3/2 ^100 keV Q-ball decay produces dark matter with 
a cosmologically interesting free-streaming length that 
can be tested in astronomical observations [62, 63]. 

This work was supported in part by DOE grant 
DE-FG03-91ER40662 and by the NASA grant 
NNX08AL48G. 



[1] A. Kusenko, Phys. Lett. B405, 108 (1997), hep- 
ph/9704273. 

[2] G. R. Dvali, A. Kusenko, and M. E. Shaposhnikov, Phys. 

Lett. B417, 99 (1998), hep-ph/9707423. 
[3] A. Kusenko and M. E. Shaposhnikov, Phys. Lett. B418, 



46 (1998), hep-ph/9709492. 
[4] M. Laine and M. E. Shaposhnikov, Nucl. Phys. B532, 

376 (1998), hep-ph/9804237. 
[5] R. Banerjee and K. Jedamzik, Phys. Lett. B484, 278 

(2000), hep-ph/0005031. 



■5 



K. Enqvist and J. McDonald, Phys. Lett. B425, 309 [34 

(1998) , hep-ph/9711514. 

K. Enqvist and J. McDonald, Nucl. Phys. B538, 321 [35 

(1999) , hep-ph/9803380. 

J. A. Frieman, G. B. Gelmini, M. Gleiser, and E. W. [36 
Kolb, Phys. Rev. Lett. 60, 2101 (1988). 
J. A. Frieman, A. V. Olinto, M. Gleiser, and C. Alcock, [37 
Phys. Rev. D40, 3241 (1989). 

K. Griest and E. W. Kolb, Phys. Rev. D40, 3231 (1989). [38 
K. Griest, E. W. Kolb, and A. Massarotti, Phys. Rev. 
D40, 3529 (1989). [39 
K. Enqvist and J. McDonald, Phys. Rev. Lett. 81, 3071 
(1998), hep-ph/9806213. [40 
K. Enqvist and J. McDonald, Phys. Lett. B440, 59 

(1998) , hep-ph/9807269. [41 
A. Kusenko, Phys. Lett. B404, 285 (1997), hep- [42 
th/9704073. [43 
A. Kusenko, Phys. Lett. B406, 26 (1997), hep- 
ph/9705361. [44 
A. Kusenko, M. E. Shaposhnikov, P. G. Tinyakov, and [45 
I. I. Tkachev, Phys. Lett. B423, 104 (1998), hep- 
ph/9801212. [46 
A. Kusenko, V. Kuzmin, M. E. Shaposhnikov, and P. G. 
Tinyakov, Phys. Rev. Lett. 80, 3185 (1998), hep- [47 
ph/9712212. 

K. Enqvist and J. McDonald, Phys. Rev. Lett. 83, 2510 [48 

(1999) , hep-ph/9811412. 

M. Axenides, S. Komineas, L. Perivolaropoulos, and [49 
M. Floratos, Phys. Rev. D61, 085006 (2000), hep- 
ph/9910388. [50 
R. Battye and P. Sutcliffe, Nucl. Phys. B590, 329 (2000), [51 
hep-th/0003252. 

R. Allahverdi, A. Mazumdar, and A. Ozpineci, Phys. [52 
Rev. D65, 125003 (2002), hep-ph/0203062. 
K. Enqvist and A. Mazumdar, Phys. Rept. 380, 99 [53 
(2003), hep-ph/0209244. 

M. Dine and A. Kusenko, Rev. Mod. Phys. 76, 1 (2004), [54 
hep-ph/0303065. 
A. Kusenko, L. Loveridge, and M. Shaposhnikov, Phys. [55 
Rev. D72, 025015 (2005), hep-ph/0405044. 
A. Kusenko, L. C. Loveridge, and M. Shaposhnikov, [56 
JCAP 0508, 011 (2005), astro-ph/0507225. 
M. Berkooz, D. J. H. Chung, and T. Volansky, Phys. [57 
Rev. Lett. 96, 031303 (2006), hep-ph/0510186. 
M. Berkooz, D. J. H. Chung, and T. Volansky, Phys. [58 
Rev. D73, 063526 (2006), hep-ph/0507218. 
A. Kusenko and A. Mazumdar, (2008), 0807.4554. [59 
M. C. Johnson and M. Kamionkowski, (2008), 0805.1748. 
S. Kasuya, M. Kawasaki, and F. Takahashi, (2008), [60 
0805.4245. 

N. Sakai and M. Sasaki, Prog. Theor. Phys. 119, 929 
(2008), 0712.1450. [61 
L. Campanelli and M. Ruggieri, Phys. Rev. D77, 043504 
(2008), 0712.3669. [62 
S. Kasuya and M. Kawasaki, Phys. Rev. D62, 023512 [63 

(2000) , hep-ph/0002285. 



M. Kawasaki, K. Konya, and F. Takahashi, Phys. Lett. 
B619, 233 (2005), hep-ph/0504105. 

S. Kasuya and F. Takahashi, JCAP 0711, 019 (2007), 
0709.2634. 

I. M. Shoemaker and A. Kusenko, Phys. Rev. D78, 
075014 (2008), 0809.1666. 

L. Campanelli and M. Ruggieri, Phys. Rev. D80, 036006 
(2009), 0904.4802. 

A. Kusenko and I. M. Shoemaker, Phys. Rev. D80, 
027701 (2009), 0905.3929. 

I. M. Shoemaker, Phys. Rev. D80, 031702 (2009), 
0907.0269. 

A. Kusenko, A. Mazumdar, and T. Multamaki, (2009), 
0902.2197. 

D. B. Kaplan, Phys. Rev. Lett. 68, 741 (1992). 

I. Affleck and M. Dine, Nucl. Phys. B249, 361 (1985). 
S. Kasuya and M. Kawasaki, Phys. Rev. D64, 123515 
(2001), hep-ph/0106119. 

F. D. Steffen, JCAP 0609, 001 (2006), hep-ph/0605306. 
R. Allahverdi and A. Mazumdar, JCAP 0610, 008 

(2006) , hep-ph/0512227. 

R. Allahverdi and A. Mazumdar, Phys. Rev. D76, 
103526 (2007), hep-ph/0603244. 

O. Seto, Phys. Rev. D73, 043509 (2006), hep- 
ph/0512071. 

S. Kasuya and M. Kawasaki, Phys. Rev. D61, 041301 
(2000), hep-ph/9909509. 

L. Roszkowski and O. Seto, Phys. Rev. Lett. 98, 161304 

(2007) , hep-ph/0608013. 

M. Ibe and R. Kitano, JHEP 08, 016 (2007), 0705.3686. 
J. L. Feng, B. T. Smith, and F. Takayama, Phys. Rev. 
Lett. 100, 021302 (2008), 0709.0297. 
A. de Gouvea, T. Moroi, and H. Murayama, Phys. Rev. 
D56, 1281 (1997), hep-ph/9701244. 

T. Moroi, H. Murayama, and M. Yamaguchi, Phys. Lett. 
B303, 289 (1993). 

T. Gherghetta, G. F. Giudice, and A. Riotto, Phys. Lett. 
B446, 28 (1999), hep-ph/9808401. 

G. F. Giudice and R. Rattazzi, Phys. Rept. 322, 419 

(1999) , hep-ph/9801271. 

M. Pospelov, Phys. Rev. Lett. 98, 231301 (2007), hep- 
ph/0605215. 

K. Jedamzik, M. Lemoine, and G. Moultaka, JCAP 

0607, 010 (2006), astro-ph/0508141. 

V. K. Narayanan, D. N. Spergel, R. Dave, and C.-P. Ma, 

(2000) , astro-ph/0005095. 

R. Barkana, Z. Haiman, and J. P. Ostriker, (2001), astro- 
ph/0102304. 

M. Viel, J. Lesgourgues, M. G. Haehnelt, S. Matarrese, 
and A. Riotto, Phys. Rev. D71, 063534 (2005), astro- 
ph/0501562. 

E. A. Baltz and H. Murayama, JHEP 05, 067 (2003), 
astro-ph/0108172. 

R. F. G. Wyse and G. Gilmore, (2007), 0708.1492. 
D. Coe and L. Moustakas, (2009), 0906.4108. 



